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We make systematic measurements of Raman anti-Stokes/Stokes (aS/S) ratios using two different 
laser excitations (514 and 633 nm) of rhodamine 6G (RH6G) on dried Ag colloids over a wide range 
of temperatures (100 to 350 K). We show that a temperature scan allows the separation of the 
contributions to the aS/S ratios from resonance effects and heating/pumping, thus decoupling the 
two main aspects of the problem. The temperature rise is found to be larger when employing the 633 
nm laser. In addition, we find evidence for mode specific vibrational pumping at higher laser power 
densities. We analyze our results in the framework of ongoing discussion on laser heating/pumping 
under SERS conditions. 
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I. INTRODUCTION 

Despite rapid progress on Surface Enhanced Raman 
Scattering (SERS)^!,,^] towards applications (S, iSfc iSl] j 
many fundamental questions remain open. An important 
issue -that has been the subject of controversy- relates to 
whether the local enhanced electromagnetic (EM) fields 
cause local heating and/or vibrational pumping. This is 
important for reasons that range from the purely practi- 
cal (damage to analytes under study) to the more funda- 
mental (vibrational pumping). In certain cases, plasmon- 
related heating effects may even be usefulQ for biomed- 
ical applications, but both the magnitude of the effect 
and how to control it need proper characterization. It is 
possible to use the aS/S ratio under SERS conditions to 
study heating/pumping effects with certain hmitations, 
as we discuss in detail below. 

The laser power dependence of the aS/S- ratio was 
thought to provide strong evidence for SERS pumping 
m , but this has been the source of controversy and alter- 
native explanations have been offered |9l llOl nVL . Con- 
tributions to the aS/S-ratio can be complex and include: 
heating, pumping, and resonance effects. There are var- 
ious hints for heating in SERS and the anomalous ratio 
(commonly measured even at low powers) was predicted 
to relate to hidden resonances in the system[lO|,lll|. We 
have recently demonstrated that it is indeed possible to 
use the aS/S ratio to map out roughly the underlying 
resonance in certain cases p^. 



Here we develop the use of aS/S-ratios further; we ex- 
pand and provide additional evidence to recent discus- 
sions in Refs. [l^ IT^ . We demonstrate how to separate 
heating from resonance effects through temperature (T) 
dependent data. We find a temperature rise at low laser 
powers (~ 0.5 mW on a 1.5 fim (diameter) spot) that 
depends on laser energy. We find evidence also suggest- 
ing the existence of mode specific vibrational pumping at 
higher laser power densities and low temperature. 



II. GENERAL CONSIDERATIONS 

For the sake of clarity we shall repeat the essential de- 
tails of the discussion presented in Refs.|0,Q| keeping 
details to a minimum. We call T^nv the temperature of 
the local environment (substrate -I- ambient surrounding 
the sample), while Troom is room temperature, Tmet the 
temperature of the metal colloids, and Tmoi the temper- 
ature of the molecule. The main possible scenarios are: 

1. No heating. The whole system is in thermal equi- 
librium, i.e. Tj-oojn — T^env — T^met — T^vnol- 

2. Global heating: system in thermal equilibrium at 



a T higher than Troom] i-e. Tei 
T 
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Local heating of substrate and analyte: when 
Tmoi = Tmet > T^nv > Troom- The aualytc is in 
thermal equilibrium with the metal, at a higher 
temperature than their environment. The source 
of heating here is mainly laser absorption in the 
metal substrate. 
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4. Local heating of analyte Tmoi ^ Tmet > T^nv ^ 
Troom- Here the molecule remains in internal ther- 
mal equilibrium at a temperature Tmoi , higher than 
the rest of the system. Strong interactions among 
vibrational modes within the molecule allows vibra- 
tional energy to thermalize and reach an equilib- 
rium through intramolecular vibrational relaxation 
or redistribution (IVR). The source of heating here 
can be direct laser absorption in the molecule (via 
electronic state transitions), or vibrational pump- 
ing (via SERS processes) followed by fast IVR. 
However, it is difficult to distinguish between these 
two possible causes. 

5. Mode specific vibrational pumping, when the an- 
alyte is not in internal thermal equilibrium. Each 
vibrational mode can be described by an effective 
mode temperature, T^^^ , which reflects the popu- 
lation of its first vibrational excited state [l^l- In 
general, IVR ensures that T^^^ is the same for all 
modes (internal thermal equilibrium at a tempera- 
ture Trnoi)- However, SERS processes affect the vi- 
brational population and can pump vibrations from 
the ground to the first excited vibrational state (for 
Stokes processes). If for a given vibrational mode, 
the pumping rate becomes comparable or higher 
than the IVR rate, this mode can then have a 
higher effective temperature than the other modes: 

^ Tfjioi ^ ^met ^ T^env ^ T^room^ this is re- 
ferred to as mode specific pumping. 

The different scenarios deserve a brief comment. If vi- 
brational pumping occurs in the presence of fast IVR, 
pumping will only be observed as an effective heating of 
the molecule (4*'* scenario) . From the experimental point 
of view, under these circumstances, it would be impos- 
sible to distinguish this situation from the 3'''' scenario. 
It is indeed very difficult to determine independently the 
temperature of the metal and the molecule, due to the 
lack of a spectroscopic feature coming explicitly from the 
metal. Even if one could prove that the 4*'' scenario was 
the correct one, it would then be extremely difficult to 
determine the real cause of heating: direct laser absorp- 
tion or vibrational pumping, or a combination of both 
effects. Consequently, the only clear-cut demonstration 
of vibrational pumping is, in our view, the last scenario; 
i.e. mode specific pumping. Heating and pumping under 
SERS conditions have previously been investigated using 
peak parameter correlations, such as width or frequency 
vs. intensity 0, Evidence was found for local heat- 
ing (3'''' or 4*'' scenario) but no evidence of mode specific 
pumping was uncovered ■ Although yielding evidence 
for the existence of local heating, it was impossible to dis- 
tinguish its origin; i.e. either within the molecule itself 
or the colloid (and consequently the molecule through 
its close contact with the metal). These previous results 
do not exclude the presence of pumping {vide supra) but 
they do not demonstrate it either. 



Most other studies of heating/pumping under SERS 
conditions 0, 0, Hol ITll | have focused on the measure- 
ments of the aS/S ratio, p. This ratio is in principle 
directly related to the vibrational populations, and there- 
fore to the effective mode temperature. However, in the 
case of SERS, it can be modified by underlying reso- 
nances associated with the wavelength dependence of the 
SERS enhancement factor, making the interpretations 
more difficult. The most general expression for the ratio 
under SERS conditions at low powers (in the absence of 
pumping) and for molecules in thermal equilibrium at a 
temperature Tmoi) is [T^ll7| : 

P = -j^ = A^'^^'S^ ex.p{-?uvv/kBT„ioi) (1) 

? — 1 

with ~ {u!l + uJv)'^/{ujl — w^)* (accounting for the 
standard wavelength dependence of Raman processes). 
Ias(S) is the intensity of the anti-Stokes (Stokes) Ra- 
man mode, ujl(v) is the frequency of the laser (Ra- 
man mode), (Jas(S) is the anti-Stokes (Stokes) scattering 
cross section of the molecule, N is the number of active 
molecules, and ^yis(s) is the SERS enhancement factor 
at the molecule position at the anti-Stokes(Stokes) fre- 
quency. The latter term includes both electromagnetic 
and chemical enhancement effects. In general the en- 
hancement can be site dependent and therefore should 
be summed over all active sites. Because anti-Stokes sig- 
nals are typically small, their measurement requires rela- 
tively long integration time and/or large analyte concen- 
trations. The results are therefore an average over several 
active sites/molecules, each with different enhancement 
factors and resonances. To account for this averaging, p 
for a given mode i can be modelled in a first approxima- 
tion as: 

p, = A^A, expi-hcuJkBT;;^^) (2) 

where Ai is the average asymmetry factor from resonance 
effects (asymmetries in cross sections) , which can change 
from mode to mode. T^'^'^ was used in this expression 
to account for the possibility of different effective mode 
temperatures. When the vibrational pumping rate is 
small compared to IVR rate, we have T^^-^ = Tmoi for 
all modes and the molecule is in thermal equilibrium at 
Tmoi- However, if vibrational pumping is not negligible, 
this expression needs to be modified to give (for weak 
pumping) [mSIIl: 

Pi = Au,Ai [e:iip{-hujJkBTmoi) + crsTrlh] (3) 

where is the vibrational relaxation time (including IVR 
and other relaxation processes), and II is the laser inten- 
sity in photons per unit time and surface area. We see 
in this latest expression that vibrational pumping would 
tend to increase the value of pi. This could alternatively 
be described by Eq. |(2Jl, using an effective mode temper- 
ature T^^^ > Tmoi: but determining T^^^ would require 
the use of Eq. ||3Jl. 
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Eqs. (0) or Q form the basis for the analysis of most 
SERS experiments about heating/pumping that use aS/S 
ratio measurements. From them, it is obvious that the 
measurement of pi only is not sufficient to extract T„ioi 
or T^^^ ^ as the AiS are unknown. One therefore needs 
an extra free parameter. In previous works, the power 
dependence of pi was for example studied H, 0, 0|. 
A linear dependence of pi with II has been invoked as 
evidence of vibrational pumping which would appear to 
follow from Eq. This is not true however and direct 
heating of the molecule has been proposed as an alterna- 
tive explanation • To understand the effect of heating 
on p, one can write Tmoi = ?o + AT, and expand Eq. Q 
considering a small temperature increase, AT ^ Tq: 

If heating effects are present, AT changes with power. 
Because heat diffusion and transfer models are linear, it 
is reasonable to assume that the temperature increase 
AT is proportional to laser power: AT = a/^. For a 
small temperature increase, the argument in the second 
exponential of Eq. is small and expanding it leads to 

p. = A^A, exp f-^] (l + ^^lA. (5) 

The linear dependence of pi (or quadratic dependence of 
I as) with incident power can therefore equally be the 
result of conventional heating effects or of pumping. It 
is also worth noting that the argument in the second ex- 
ponential of Eq. can be comparable to 1, even for 
a small AT. This can be the case for example for high 
energy modes, where TiWi/ZcsTo can be « 10. The power 
dependence of pi as a result of heating effects would then 
be exponential. Such power dependence, linear for low 
energy modes, and exponential for high energy modes, 
has already been observed [31 • This shows that heating 
effects can be important and that a linear power depen- 
dence of Pi is not in itself a proof of vibrational pumping. 

Instead of using the laser power as an additional pa- 
rameter, we propose to study the aS/S ratios of each 
mode as a function of substrate temperature Tq. As we 
shall show, by investigating how these ratios change as 
a function of T over a wide range, we can separate the 
contributions due to heating and resonance effects. In 
addition, this approach enables us to detect a departure 
from internal thermal equilibrium of the molecule, i.e. 
vibrational pumping. 

III. EXPERIMENTAL DETAILS 

Silver colloids were produced as described in Ref.fTsf. 
SERS samples were prepared by mixing equal amounts 
of colloidal and 25 mM KCl solutions together, to which 
1 pL of M Rhodamine 6G (RH6G) solution was 

added, resulting in a 10~^ M RII6G concentration. The 



sample was then dried onto a Si substrate. The tempera- 
ture calibration in terms of aS/S-data of non-SERS active 
signals (paracetamol) is described elsewhere fl^. Mea- 
surements were carried out with a Renishaw 2000 CCD 
spectrometer and an Olympus BH-2 microscope with a 
Linkham temperature stage, using both 514 nm Ar+ ion 
and 633 nm HeNe lasers with a x50 long working dis- 
tance objective (^^1.5 /im beam diameter at the focal 
point). Measurements with the 514 nm laser were made 
at 0.5 mW whereas for the 633nm laser were performed 
at 0.5 and 5 mW, respectively, in the range between 100 
and 350 K in steps of 10 K and integration times varying 
from 1 to 30 sec. Five or more measurements were made 
to gain an average over the sample. Peaks were then 
analyzed using standard Voigt functions with subtracted 
backgrounds. 

The variation in the aS/S ratio for each mode was fitted 
to Eq. IP, with T^^^ = Ts, -t- AT,., where AT, is the in- 
crease in effective mode temperature with respect to the 
T of the substrate which is measured by the aS/S-ratio 
of the Si substrate^^. The parameters Ai in Eq. Q 
will depend on the particular mode under consideration, 
accounting for the asymmetry in 175 and Oas produced 
by resonances. Four types of fits were performed: 

• AT-only, where Ai is set to 1 (no resonance effect) 
and ATi arc fitting parameters; 

• A-only, where AT^ = (no local heating, no pump- 
ing) and Ai are fitting parameters; 

• A+AT-shared, where AT^ = AT constant for all 
modes and Ai are fitting parameters. This corre- 
sponds to the case of internal thermal equilibrium 
for the molecule; 

• A4-AT, where Ai and AT; are independent fitting 
parameters. This allows the modes to have different 
effective mode temperatures. 

Figure ^ illustrates the predictions of the model and 
the influence of A and AT on the value of the aS/S-ratio 
as a function of T. Figure^^a) shows the expected value 
of the ratio for different A's, whilst keeping AT = 0. 
When ^ < 1 (A > 1) the ratio is suppressed (increased). 
In Fig^b) A is held constant whilst AT is varied: the in- 
fluence of AT is strongest at low T's where it is a greater 
fraction of the initial T, resulting in a distinctive 'flar- 
ing' of the ratio at low T's and a flattening of the overall 
curve shape. 

The flrst two types of fits are only used to test whether 
a simple model with heating only and no resonances, 
or resonances only and no heating, is sufficient to ex- 
plain the data. For all fits, we assume that Ai and AT^ 
are constant with temperature. Ai, which is the result 
of predominantly electromagnetic resonances, is not be- 
lieved to vary with temperature. For molecules in ther- 
mal equilibrium AT; — AT is the same for all modes and 
represent the temperature increase (heating). Many fac- 
tors contribute to this increase, for example the thermal 
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A + AT shared 


(A + AT) 


A only 


AT only 


Mode (cm"^) 


A 


A 


AT 


A 


AT 


Low power 633 nm laser 


610 


1.30 ± 0.05 (AT = 25K) 


1.50 ± 0.05 


17 ± 5 


2.00 ± 0.20 


40 ± 5 


780 


1.35 ± 0.05 (AT = 25K) 


1.65 ± 0.10 


16 ± 5 


2.35 ± 0.20 


35 ± 5 


1310 


0.85 ± 0.10 (AT = 25K) 


0.50 ± 0.10 


55 ± 10 


1.60 ± 0.30 


20 ± 5 


1360 


0.95 ± 0.10 (AT = 25K) 


0.80 ± 0.10 


35 ± 5 


1.90 ± 0.30 


25 ± 5 


1510 


1.05 ± 0.10 (AT = 25K) 


0.70 ± 0.10 


40 ± 6 


2.20 ± 0.40 


30 ± 10 


High power 633 nm laser 


610 


0.45 ± 0.05 (AT = 96K) 


1.05 ± 0.10 


40 ± 5 


3.1 ± 0.5 


40 ± 5 


780 


0.40 ± 0.04 (AT = 96K) 


0.70 ± 0.10 


60 ± 15 


4.4 ± 1 


50 ± 5 


1310 


0.50 ± 0.10 (AT = 96K) 


0.15 ± 0.05 


160 ± 35 


12 ± 3 


75 ± 10 


1360 


0.40 ± 0.10 (AT = 96K) 


0.25 ± 0.10 


120 ± 35 


14 ± 3 


70 ± 10 


1510 


0.50 ± 0.10 (AT = 96K) 


0.08 ± 0.03 


200 ± 40 


20 ± 6 


80 ± 10 


Low power 514 nm laser 


610 




0.45 ± 0.05 


3 ± 2 


0.50 ± 0.03 


-20 ± 4 



TABLE L Fitting results for the 514 and 633 nm lasers using the models described in the text. Results are shown for both 
high and low power measurements for the 633 nm laser whilst only low power measurements were possible with the 514 nm 
laser. Low power means ~0.5 mW at the focal point whilst the high power measurement was an order of magnitude larger. 
All measurements were made on large colloidal clusters. The j4-|-Ar-shared fitting is not possible for the 514 nm laser as only 
one mode was measurable. The (A+AT-shared) and (A+AT) models provide the best representation of the data at low and 
high power respectively. 



properties of the SERS substrate and of the various in- 
terfaces. These factors could possibly vary with temper- 
ature, although large variations typically occur mainly 
at even lower temperatures (T < 100 K). The assump- 
tion AT =constant is therefore reasonable in a first ap- 
proximation. The (A-|-AT-shared) fit then enables us to 
separate the resonance effects (Ai) from the heating ef- 
fect (AT). If this fit fails, then either AT ^constant is 
wrong, or the molecule is not in internal thermal equi- 
librium (i.e. pumping occurs). In the latter case, the 
effective mode temperature Tq + ATi is determined by 
the pumping rate, via Eq. (Q. The fits with this model 
cannot therefore be valid, but the (A+AT) fit can still be 
useful in identifying which modes present an anomalous 
effective mode temperature and are therefore most likely 
to be pumped. 

If the (74-l-AT-shared) fit fails, the model above can- 
not distinguish directly between pumping (AT^ different 
for different modes) and temperature dependent heating 
(ATi = AT, but AT varies with T). In order to prove 
that different modes cannot be accounted for by a single 
T, we propose to define the following function (for two 
modes a and b): 



F 



fcBln(pa) kBlnipb) 



huja 



hujh 



(6) 



Then, describing each mode by an effective mode tem- 
perature, we have from Eq. 



F = C 



1 



T 



(7) 



where C is a constant accounting for the different A- 
factors. The advantage of this function is as follows: 
if both modes share the same temperature, this func- 
tion is constant as a function of T. The function F can 
be used to monitor relative changes between the effec- 
tive T's of two modes as a function of external varying 
conditions (here T, but it could also be used to prove 
non-equilibrium pumping in a power dependence experi- 
ment). 



IV. RESULTS AND DISCUSSION 
A. Low power 

Figure|21shows examples aS/S ratios as a function of T 
for the 610 cm~^ mode of RH6G and 633 nm excitation 
(top curve). The theoretically predicted behavior (A=l, 
AT = 0) and the fits to the experimental data using the 
[A+AT], A-only, and AT-only models, are shown. The 
experimental data are well approximated by the (A-f-AT- 
shared) fit for all modes investigated. The fitting to all 
the modes, using the four fitting procedures presented 
above, is summarized in tabled Figure |21 also shows the 
equivalent data for the 610 cm~^ mode under the 514 nm 
laser excitation (bottom curve). The aS-peaks associated 
with the higher modes were too small to analyze at these 
laser powers. The (A -I- AT) fitting provides the best fit 
to the experimental data although AT is small and con- 
sequently the A only model also fits the data reasonably 
well. The fitting results for the 610 cm~^ mode taken 
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using the 514 nm laser are also summarized in tabled 

The low power measurements tell us that the molecule 
is in thermal equilibrium with a heating of 25K when us- 
ing 633nm laser and 3K when using the same low power 
514nm laser. Useful insight into the colloidal heating 
problem can be gained from considering a simplified 1-D 
model. Consider a single colloid where the laser power 
is assumed to dissipate through the top surface with 
no lateral heat flow. The T dependence of the ther- 
mal conductivity of Ag is small down to 40 K and we 
can use the room temperature value. The internal ther- 
mal time constant of the colloid (the time for the col- 
loid to come to a uniform temperature), is given by 
Tc = w 1 X 10^^^ (s), where Kc is the thermal con- 
ductance of a single Ag particle (of size ~ (20 x 10^^ 
m)'^) and Caq is the phononic contribution to the spe- 
cific heat with Debye temperature, QAg given by: Caq = 

CAg {{DAg/MWAg X lO'^) X (P) » 6.1 X IQ~^^{J/K). 

In this latter expression WAg is the atomic weight of 
Ag, DAg its density, d the length of the colloid and 
CAg ~ 7.8 X 10^ {J / M ol .K) . We ignore the electronic con- 
tribution to CAg as it is much smaller than the phonon 
contribution, except at liquid He temperatures. 

Assuming an incident power 1 mW is focussed into a 1 
jjLw? area with a reflectivity of 80%, the number of pho- 
tons absorbed per second by the colloid is « 2.6x lO^^s"^. 
Note that when dealing with small objects, Tc is compa- 
rable to the arrival rate of the photons. The dynamics 
of the system is therefore rather interesting. The ab- 
sorption of a single photon increases the colloid's upper 
surface temperature by « 5.2 x 10^"^ K. By modelling 
the acoustic mismatch between Ag and Si 19], we can 
estimate the expected increase in T. Estimating val- 
ues from experimental data for interfaces between Al- 
AI2O3 and Pb-diamond, we find temperature increases 
of 2 and 10 K, respectively, which are the of the right 
order of magnitude for our low power observations using 
the 633nm laser. This is remarkable when we consider 
that we have ignored any potential contribution origi- 
nating from SERS hot-spots and that the thermal model 
is extremely simplified. This model supports the sce- 
nario that at low powers the molecules could be heating 
through direct interaction with the colloids. 

B. High power 

Figure 13 shows p as a function of T for the 610 and 
1510 cm""'^ modes for high power (5 mW) 633 nm exci- 
tation, while all other parameters were identical to the 
measurements at low power. We show the resulting fits of 
[A+AT] for both modes, as well as secondary fits swap- 
ping the temperatures and the theoretically expected p's. 
The high power ratios for all the modes were also fit- 
ted using the (A-|-AT-shared) model as in the low power 
data. This resulted in a AT = 96K whereas the low 
power results gave 25 K. Although this is clearly larger, 
as would be expected, the quality of the fits were badly 



compromised for the high power data, as measured by a 
X^-test (not shown in the table) . This indicates that this 
model is no longer valid in this regime. The results of the 
fitting to the high power, 633 nm data for all the modes 
is also summarized in table The 610 cm^^ ratio as a 
function of T at high power is similar to that obtained at 
low powers; when plotted on the same graph it is almost 
impossible to separate them by eye suggesting that the 
mode temperature is not too dissimilar, despite an order 
of magnitude difference in the incident power. Indeed, 
the result of the {A+/S.T) fitting of the 610 cm^^ mode 
in both cases gives AT — 17±5 K for low power and AT 
= 40 ± 5 K for high power, which is difficult to visualize 
on a log-scale. For the 1510 cm~^ mode, although there 
is a great deal of scatter, there is a noticeable difference 
in the trend of the data when compared with the low 
power measurement. This behavior is very similar to the 
'fiaring' of the theoretically determined p's for high AT 
as shown in figure ^ This is confirmed with the fitting of 
the data to the (A-I-AT) model which results in AT1510 
= 200 ± 40 K; i.e. much larger than the Tgio = 40 ± 5 K 
from the best fit to the 610 cm~^ mode at high power. 
FigureOalso shows fits for the 610 and 1510 cm~^ modes 
with ATeio =200K and AT1510 = 40, both with the best 
fit value of A for that mode. It is obvious from these fits 
that the 610 cm"-'^ mode is much cooler than the 1510 
cm~^. All these results suggest that vibrational pump- 
ing in the 1510 cm^^ mode could be the explanation for 
the failure of the (A+AT-shared) fits. 

In order to confirm this hypothesis, we plot in Fig. 01 
the function T, discussed in the previous section, which 
compares the effective T's of (a) the 780 and 610 cm~^ 
modes and (b) the 1510 and 610 cm~^ modes in (b). The 
monotonic departure from a constant observed in (b) is 
an additional demonstration of non thermal equilibrium 
between these two modes and therefore in the molecule. 
This confirms that the failure of the (v4-f AT-shared) fits 
is a consequence of vibrational pumping, and not of a 
temperature dependent heating (same AT for all modes, 
but varying with temperature). In fact, the data in Fig. 
21 show not only evidence for mode specific pumping but 
also highlights the limitations of the A-I-AT model. The 
kink in T as a function of T means that the 1510, 780, 
and 610 cm~^ modes share a single T above ^ 200 K, 
whilst below that the 1510 cm"-'^ shows some mode spe- 
cific pumping. This is only taken into account approxi- 
mately in the A~\-AT model, and Eq. lO needs to be used 
to account for pumping properly. Fits using Eqs. Q and 
© are also shown in Fig. 0] and are in good agreement 
with the data. Moreover, we extract from the fit for the 
1510 cm^^ mode a value of crsTriL = 5 x 10"^. We es- 
timate II ~ 10^^ photons/cm^ and typical values for 
are in the range 1-10 ps. This leads to as in the range 
5x 10~^^— 5x 10^^^, compatible with previous estimate of 
the SERS cross section. The study of the T-dependence 
of the function F therefore provide additional evidence 
that the behavior of the 1510 cm~^ below 200 K is a mani- 
festation of vibrational pumping in this mode. This effect 
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is more prominent at low T and for higher energy modes 
because the equihbrium population of the excited vibra- 
tional state is then very small and the relative change 
in intensity is larger. Following this argument, further 
evidence for pumping could be obtained from measure- 
ments at even lower temperatures (T < lOOK) than the 
ones accessed in here. We note that the original paper 
on SERS pumpingalso presented some evidence for mode 
dependent effects Our results certainly contribute to 
this body of work and suggest that several mechanisms 
(vibrational pumping, underlying resonances, and colloid 
or SERS substrate heating) operate simultaneously and 
contribute together to the aS/S ratio anomalies. 

In summary, we detail the key results presented in ta- 
ble as follows: (i) The (A-|-AT-shared) model provides 
a satisfying agreement at low powers. The molecule is 
then in thermal equilibrium with a heating of '--^ 25 K 
with respect to the substrate for the 633 nm laser, but 
only a 3 K rise for 514 nm excitation, (m) Models with 
only A or AT alone as fitting parameters systematically 
produce worse overall agreement with the data or even 
unphysical values, like negative heating (see for example 
the 610 cm^^ mode at 514 nm excitation in tableP). {Hi) 
For high powers using the 633 nm laser, an overall shared 
AT among the modes is not possible and one is forced 
to a better quality fit by allowing different effective tem- 



peratures for different modes. Using the function F, one 
can then demonstrate that this is the result of vibrational 
pumping in the 1510 cm~^ mode. 



V. CONCLUSION 

In closing, we have observed a laser frequency depen- 
dent heating of the order of 25 K at low laser power under 
SERS conditions with 633 nm excitation. We have shown 
that it is possible to decouple in dry samples the contri- 
butions of the thermal population and resonances in the 
aS/S ratios by performing a T-scan, thus overcoming the 
basic uncertainty in the different relative contributions 
which exist in studies at a single T. Most importantly we 
have used this technique to provide evidence for at least 
two possible scenarios: molecule in thermal equilibrium 
at low powers and mode specific pumping at high powers. 
Experiments at lower temperatures are in progress and 
will be reported elsewhere. 
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FIG. 1: Theoretical plots showing the efltoct on the anti- 
Stokes/Stokes ratio of various values of A (left) and AT 
(right). Differing values of A result in nested anti- 
Stokes/Stokes ratios over the full temperature range. In the 
case of AT, the effect on the anti- Stokes/ Stokes ratio varies 
across the full temperature range with the largest effect being 
made at lower temperatures where AT represents a greater 
percentage of the initial temperature. This results in a 'flar- 
ing out' of the ratio at lower temperatures which is noticeably 
larger for greater values of AT. 
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FIG. 2: Anti-Stokcs/Stokcs ratios as a function of tempera- 
ture for the 610 cm^'^ SERS active modes of RH6G for the 
633 nm and 514 nm lasers at low power (0.5 mW). There is a 
clear different in the resonance contribution to the anomalous 
ratios between these lasers as the 633 nm ratio is increased 
above the theoretical value whilst the 514 nm ratio is sup- 
pressed below the theoretical. The solid lines represent the 
fit to the experimental data using the (A+AT) model whilst 
the dashed and dash dot represent the fits obtained from us- 
ing only the AT or A parameter respectively. The dotted line 
is the theoretical ratio values. See text for details. 
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FIG. 3: Anti-Stokcs/Stokcs ratio as a function of temperature 
for the 610 and 1510 cm^^ SERS active modes of RH6G taken 
with the 633 nm laser at high power (5 mW). The theoreti- 
cally expected ratio is shown as the dotted line. Two fits with 
the same value of A are shown for each mode. ATeio = 40K 
and ATisio = 200K represents the best fit values for these 
modes. The secondary fits set ATeio = 200K and AT1510 = 
40K. 
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FIG. 4: F function (Eq. vs. T comparing the effective 
temperatures of two different pairs of modes: (a) 780 and 610 
cm~^ and (b) 1510 and 610 cm~^. A constant value of F 
vs. T in (a) means that these two modes share a common 
temperature. The comparison in (b) however points towards 
mode selective pumping of the 1510 cm~^ mode with respect 
to the 610 cm~^. Solid lines are fits to the data. See the text 
for further details. 



